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ABSTRACT 

Long baseline oscillation experiments may well emerge as test beds for neutrino inter- 
actions as are present in R-parity violating supersymmetry We show that flavour diagonal 
(FDNC) and flavour changing (FCNC) neutral currents arising therefrom prominently im- 
pact a neutrino /3-beam experiment with the source at CERN and the detector at the pro- 
posed India-based Neutrino Observatory. These interactions may preclude any improvement 
of the present limit on #13 and cloud the hierarchy determination unless the upper bounds on 
couplings, particularly A', become significantly tighter. If Ifi interactions are independently 
established then from the event rate a lower bound on # 13 may be set. We show that there 
is scope to see a clear signal of non-standard FCNC and FDNC interactions, particularly in 
the inverted hierarchy scenario and also sometimes for the normal hierarchy. In favourable 
cases, it may be possible to set lower and upper bounds on A' couplings. FCNC and FDNC 
interactions due to A type Ijt couplings are unimportant. 
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I Introduction 



Neutrino physics, a bit player on the physics stage in yesteryears, has now donned a cen- 
tral role. Various experiments on atmospheric pQ, solar [2], reactor [Hj, and long baseline 
neutrinos [I] not only indicate oscillation but also pin down most neutrino mass and mixing 
parameters - two mass-square differences and two mixing angles. The best-fit jS] values with 
3cr error from atmospheric neutrinos are |Am| 3 | ~ 2. 12 jr^si x e ^ 2 , #23 — 45.0°+g° 3 p 
and from solar neutrinos Am 2 2 ~ 7.91^ x 10~ 5 eV 2 , 9 12 ~ 33.21°^o. Here Am 2 = 
m'j — mf. The sign of Arn| 3 is yet unknown and the neutrino mass spectrum will be referred 
here as normal (inverted) hierarchical if it is positive (negative). Using reactor antineutrinos 
[UlEl El; an upper bound has been set on the third mixing angle at 3<r level as sin 2 # 13 < 0.044 
resulting in 6* 13 < 12.1°. The phase 5 in the neutrino mixing matrix is not known. 

Research in this area is now poised to move into the precision regime. Can we use 
upcoming neutrino experiments to probe non-standard interactions like ^? supersymmetry? 
If present, can they become spoilers in attempts to further sharpen the neutrino properties? 
We attempt to address these issues with a specific experiment as our laboratory. 

Long baseline neutrino oscillation experiments using neutrino factories [7j and /3-beams 
jHl El [Jill HI] hold promise of refining our knowledge of # 13 , 5, and the sign of Am^. A 
possible experiment in this category would use neutrinos from a /3-beam from CERN to the 
proposed India-based Neutrino Observatory [12J (INO), a baseline of ~ 7152 Km. This is 
the set-up which we consider here. For such an experiment the /5-beam is required to be 
boosted to high 7. 

Interaction of neutrinos with matter affect long baseline experiments and this becomes 
more prominent at higher values of #13. Various authors ^Hl have considered this effect for 
atmospheric neutrinos. 

Apart from the electroweak effects, there may well be non-standard interactions leading 
to flavour changing as well as flavour diagonal neutral currents (FCNC and FDNC). Here we 
have in mind interactions with quarks and leptons involving an initial and a final neutrino. If 
there is no change in the neutrino flavour - as, for example, in Z° exchange - this is classified 
as an FDNC process, while it would be FCNC otherwise. R-parity violating supersymmetric 
models (RPVSM) [Tl|, which have such interactions already built-in 1 , will be the main focus 
of our work. Very recently a model in which couplings associated with FCNC and FDNC can 
be quite a bit higher than permitted in RPVSM has also been considered |15| I16|. Naturally, 
here the matter effect will be further enhanced. However, as RPVSM is a well-studied, 
renormalizable model which can satisfy all phenomenological constraints currently available, 
we shall restrict our main analysis only to it and shall make qualitative remarks about the 
other model, for which our results can be easily extended. 

Consequences of FCNC and FDNC for solar and atmospheric neutrinos ^3 EI], and 
neutrino factory experiments ^H] have been looked into. Our focus is on /9-beam experiments, 
particularly over a long distance (7152 Km) baseline. Our analysis encompasses both normal 
and inverted hierarchies and we also incorporate all relevant trilinear R-parity violating 
couplings leading to FCNC and FDNC. Huber et al [20] have a somewhat similar analysis 
using neutrino beams obtained from muon decays. 

The very long baseline from CERN to INO will capture a significant matter effect and 
offers a scope to signal non-standard interactions. We examine whether the presence of Tfc 

1 e.g. through squark (A'-type couplings) or slepton (A-type couplings) exchange. 
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interactions will come in the way of constraining the mixing angle #13 or unraveling the neu- 
trino mass hierarchy. The possiblity to obtain bounds on some R-parity violating couplings 
is also probed. 

II /3-beams 

A beta-beam jH], is a pure, intense, collimated beam of electron neutrinos or their an- 
tiparticles produced via the beta decay of completely ionized, accelerated radioactive ions 
circulating in a storage ring [21]. It has been proposed to produce v e beams through the 
decay of highly accelerated 18 Ne ions and v e from 6 He [21] ■ More recently, 8 B and 8 Li [22] 
with much larger end-point energy have been suggested as alternate sources since these ions 
can yield higher energy v e and v e respectively, with lower values of the Lorentz boost 7. It 
may be possible to have both beams in the same ring, an arrangement which will result in a 
v e as well as a v e beam pointing towards a distant target. In such a set-up the ratio between 
the boost factors of the two ions is fixed by the e/m ratios of the ions. Here, we will present 
our results with the 8 B ion (Q = 13.92 MeV and ti/ 2 = 0.77s) taking 7 = 350. As we show, 
7 ~ 350 may permit a distinction between matter effects due to Standard Model interactions 
and those from R-parity violating supersymmetry (SUSY). Details of the neutrino flux from 
such a /3-beam can be found in [TT] . 

Using the CERN-SPS at its maximum power, it will be possible to access 7 ~ 250 2 |2"3"] . 
For a medium 7 ~ 500, a refurbished SPS with super-conducting magnets or an acceleration 
technique utilizing the LHC [23 I2H EH] would be required. In the low-7 configuration, 
1.1 x 10 18 useful decays per year can be obtained with 18 Ne ions [2H1 [21] • We have used 
this same luminosity for 8 B and higher 7 [2"S| . This issue is being examined in an on-going 
dedicated machine study at CERN. 

III R-parity violating Supersymmetry 

In supersymmetric theories [T3|, gauge invariance does not imply baryon number (B) and 
lepton number (L) conservation and, in general, R-parity (defined as R = (— \ s j 3B + L + 2S where 
S is the spin) is violated. To maintain consistency with non-observation of fast proton decay 
etc, in the R-parity violating Minimal Supersymmetric Standard Model (imposing baryon 
number conservation) one may consider the following superpotential 

Wjz = / ] -j^ijkULjEl + X'^LiQjDl + [liLiEu, (1) 

i,j,k 

(suppressing colour and SU{2) indices) where i, j, k are generation indices. Here Lj and Qi 
are S77(2)-doublet lepton and quark superfields respectively; E iy Di denote the right-handed 
SU (2)-singlet charged lepton and down-type quark superfields respectively; H u is the Higgs 
superfield which gives masses to up-type quarks. Particularly, \j k is antisymmetric under 
the interchange of the first two generation indices. We assume that the bilinear terms have 
been rotated away with appropriate redefinition of superfields and focus on the two trilinear 
L-violating terms with A and A' couplings. Expanding those in standard four-component 

2 7 = 250 yields too few events in this experiment for the extraction of interesting physics. 
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Dirac notation, the quark-neutrino interaction lagrangian can be written as: 



c y = \' ijk [ 4 d* 4 + (4)*(4) c 4l + h.c. (2) 

Above, the sfermion fields are characterized by the tilde sign. The charged lepton interacts 
with the neutrino via 

£a = l^k[ei4ul + (e k R r(uirei-(t^j)]+h.c, (3) 

In what follows, we shall consider these couplings as real but will entertain both positive and 
negative values. The interactions of neutrinos with electrons and rf-quarks in matter induce 
transitions (i) v i + d — * Vj + d and (ii) v { + e — > Vj + e. (i) is possible through A' couplings 
via squark exchange for all i,j and through Z exchange for i—j while (ii) can proceed via 
W and Z exchange for i — j, as well as through A couplings via slepton exchange for all i,j. 



IV Neutrino oscillation probabilities & matter effect 

In a three neutrino framework, the neutrino flavour states |i/ a ), a = e,/j,,T, are related to 
the neutrino mass eigenstates i = 1,2,3, with masses mf 

Wa) = ^2U ai \ui) , (4) 

i 

where U is a 3 x 3 unitary matrix which can be expressed as: 

U = V 23 V 13 V 12 , (5) 



where 

Cl3 

r 2 , = ! o .s 2 , j : r l:! = | o :i o | : r L2 

s 13 e i<5 





/ C12 S12 

-Sl2 C12 | . (6) 

\ 1 



Qj = cos Oij, Sij = sin and 5 denotes the CP violating (Dirac) phase. There may be a 
diagonal phase matrix on the right containing two more Majorana phases. These are not 
considered below. Apart from some qualitative remarks, we present our result considering 
5 = corresponding to the CP conserving case. In the mass basis of neutrinos 

M 2 = diag^, m \, m\) = CM; M,J\ (7) 

where M v is the neutrino mass matrix in the flavour basis and mi, m,2, and m 3 correspond to 
masses of three neutrinos in the ascending order of magnitudes respectively for the normal 
hierarchy. If m 2 > m 3 , we have an inverted hierarchy. 
The neutrino flavour eigenstates evolve in time as: 





(8) 



where 

H = Exl 3x3 + u(—jU^ + R. (9) 

Here E is the neutrino energy while I3X3 is the identity matrix. R is a 3 x 3 matrix reflecting 
the matter effect, absent for propagation in vacuum. 



Specifically, 



Ri, = Rij(SM) + Rij(X') + RijiX). (10) 



Rij(SM) = y/2G F n e 5ij{i,j = 1) + (H) 

a/2 



(A' A' A' A' \ 

, im 2 (d m ) d Am 2 (d m ) V 



(12) 



\Am 2 {d m ) 4m 2 (d. 

Rij(\) = y ^en e + y ^^n e , (is) 

where Gp is the Fermi constant, n e , n n , and rid, respectively, are the electron, neutron, 
and down-quark densities in earth matter. Note that R is a symmetric matrix and also 
that antineutrinos will have an overall opposite sign for R^j. Assuming earth matter to be 
isoscalar, n e = n p = n n and = 3n e . The current bounds on the $ couplings (Ej imply 
that the A' induced contributions to i?n, R± 2 and Ri 3 are several orders less than \f2Gpn e . 
We neglect those terms in our analysis. The upper bounds on all couplings in Rij(X) are also 
very tight [H] in comparison to \f2Gpne and their effect will be discussed later. So, first we 
consider, in addition to the Standard Model contribution, only 

-R23 = -R32 = - — d ~ {X' 2m \X' Zml + X' 21m X' 31m ) , 
Am 2 (d m ) 

R22 = - — d ~ {X' 2 ml + Kim) 1 -^33 = ~ — " J {X' 3m i + X' 3 l m ) , (14) 
4771 yCL m ) 4771 \CLm ) 

which are comparable to V2Gpn e . One can see from eq. ((H)) that R 23 ^ implies both R 22 
and R 33 are non-zero 3 . 

The current bounds on the relevant couplings are as follows [Hj: 

|A' 221 , A' 231 | < 0.18; |A' 21 J < 0.06; |A' 331 | < 0.58; |A' 321 | < 0.52; |A 31 J < 0.12, (15) 

for down squark mass mj = 100 GeV. The chosen limits on A' 21m and A' 31m do not conflict 
with the ratio R T7T = T(r — > 7iu T )/T(7T — > fiu^) |14j . However, the recently published 
BELLE bound on the mode r — > ji-n 29j tightly constrains precisely those products of the 
A' couplings which enter in R 23 = R 32 in eq. (JHJ). It has been shown that |A' 21m A 31m | and 
|A 2ml A 3ml | both must be < 1.8 x 10~ 3 ( 100 ?T q c , v ) 2 [30J. This effectively makes R 23 negligible 
for our purposes. 

In general, it is cumbersome to write an analytical form of the probability of neutrino 
oscillation in the three-flavour scenario with matter effects. However, under certain reason- 
able approximations it is somewhat tractable. Firstly, a constant matrix can be extracted 
from M 2 in eqs. (JZj) and (|9~|) . Also for the energies and baselines under consideration, 

3 However, in other models JS] this may not be the case. 
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Am\ 2 L/ E « 1. Under this approximation, the V12 part of U drops out from eq. 0. With 
these modifications, the effective mass squared matrix is: 



2E V y/2 



(16) 



In the special case where R22 = R33, if one uses the best-fit value of the vacuum mixing 
angle 823 = vr/4 then the neutrino mass squared eigenvalues are: 



2E 



Roo — R 



23- 



IE 



^^13 

IE 



+ R u + R 22 + R23TA 



(17) 



where 



A 



Am 2 3 
2E 



+ (-R n + R22 + R23? - 2 



Am 



13 



2E 



cos 2^i3 {Ru - R22 - R23) 



1/2 



(18) 



The matter induced neutrino mixing matrix is given by 



IT 



US 




(19) 



Here 



^1, 



2 (-R n + R 22 + R 2 3 + ^§ 1 cos 20 13 ± A ) 



1 -1/2 



Am i3 
2E 



+ 2 



sin 2 29 



13 



(20) 



where N% {N3) corresponds to +(— ) sign in the above expression. Neglecting the CP phase 
in the standard parametrization of U m , one may write = sin 6*^ and = sin#^; cos#^. 
From eq. (|19j) it follows that 6^ = 023, the vacuum mixing angle. 6^, on the other hand, 
changes from its vacuum value and it is 7r/4 for 



Rn — R 



22 



R 



23 



Am?o 

— cos 20 
2E 



13- 



(21) 



In the absence of non-standard interactions, R 



22 



R 



23 



R 



33 



and R\\ = V2Gpn e , this 
is the well-known condition for matter induced maximal mixing. Since in eq. ()19|) = 0, 
in the v e to oscillation probability the terms involving (Mf — M 2 ) and (M| — M|) will 
not survive and we get: 



4 (U? 3 f {U™f sin 2 (1.27 A L) 



(22) 



where E, Am 2 3 and L are expressed in GeV, eV 2 , and Km, respectively. This expression is 
also valid for antineutrinos. Using eqs. (|17j) and (|19|) one can easily obtain the oscillation 
probabilities for other channels. 

We use the above analytical formulation as a cross-check on our numerical results. For 
example, Fig. ^ which shows the variation of P v ^ u as a function of the energy, is obtained 
using the full matter-induced three-flavour neutrino propagation including non-standard 
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Figure 1: P v ^ v for the normal and inverted mass hierarchies. SM corresponds to only standard 
electroweak interactions. The values of A' are given in parantheses. m can take any value, n = 2 or 3. 



interactions. The range of energy is chosen in line with the discussions in the rest of the 
paper. The probability falls with decreasing 13 and, for illustration, we have chosen a value 
in the middle of its permitted range. The purpose of Fig. Q is twofold: (a) to show how 
the distinguishability between the normal and inverted hierarchies may get blurred by the 
RPVSM interactions, and (b) how irrespective of the hierarchy chosen by Nature the results 
may be completely altered by the presence of these interactions. Each panel of Fig. ^ has 
three curves: the solid line (only electroweak interactions) , dot-dashed line (in addition, .R33 
gets a non-zero RPVSM contribution), and dashed line (R22 — R33 are nonzero, in addition 
to the electroweak contribution). Only in the last case is the analytical formula we have 
presented above applicable. We find excellent agreement. Two aspects of the results are 
worth pointing out. 

First, in the absence of non-Standard interactions, for an inverted hierarchy the resonance 
condition eq. (J21j) is not satisfied and the oscillation probability is negligible (right panel 
solid line). This could be altered prominently by the RPVSM interactions (dot-dashed curve) 
so that the distinguishability between the two hierarchy scenarios may well get marred by 
R SUSY. 

Secondly, for the normal hierarchy, it is seen that the peak in the probability may shift to a 
different energy in the presence of the RPVSM interactions. This is because the condition for 
maximal mixing in eq. (|21j) is affected by the R interactions. For the inverted hierarchy, the 
oscillation probability is considerably enhanced for some energies. Thus, physics expectations 
for both hierarchies will get affected by RPVSM. 

In the following section we dwell on the full impact of this physics on a long baseline 
/3-beam experiment. 



V Results 

We consider a long baseline experiment with a v e /3-beam source. /3-beams producing v e 
are also very much under consideration. Broadly speaking, the results obtained for a v e 
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beam with a normal (inverted) hierarchy are similar to that with a v e beam for an inverted 
(normal) hierarchy but details do differ. 

The average energy of the v e beam in the lab frame is (E) = 2jE cms , where E cms is 
the mean center-of-mass neutrino energy. With 7 = 350 and Q = 13.92 MeV for an 8 B 
source, (E) ~ 5 GeV. The proposed ICAL detector at INO jT2j consists of magnetized iron 
slabs with glass resistive plate chambers as interleaved active detector elements. We present 
results for a 50 Kt iron detector with energy threshold 2 GeV. As signature of v e — ► 
oscillation, prompt muons will appear 4 . Their track reconstruction will give the direction 
and energy of the incoming neutrino. ICAL has good charge identification efficiency (~ 95%) 
and a good energy resolution ~ 10% above 2 GeV. Details about the detector and neutrino 
nucleon cross sections may be found in [TT]. For the cross section we include contributions 
from quasi-elastic, single pion, and deep inelastic channels. For our chosen high threshold (2 
GeV), the contribution from the deep inelastic channel is relatively large. For the CERN- 
INO baseline, the averaged matter density is 4.21 g cm -3 . We use best-fit values of vacuum 
neutrino mixing parameters as mentioned in the Introduction. All the presented results are 
based on a five-year ICAL data sample 5 . 

At the production and detection levels, FCNC and FDNC effects can change the spectrum 
and detection cross sections by a small (< 0.1%) amount but this would not alter the 
conclusions. At the source and detector, they may also mimic the oscillation signal itself, 
but these effects are tiny 6 (~ 0(1O~ 14 )). Here we discuss how FCNC and FDNC may 
significantly modify the propagation of neutrinos through matter over large distances. 

V.l Extraction of #13 and determination of hierarchy 

If neutrinos have only Standard Model interactions then the expected number of muon events 
is fixed 7 for a particular value of #13 with either normal or inverted hierarchy as may be seen 
from the solid lines in Fig. |2| The vast difference for the alternate hierarchies picks out 
such long baseline experiments as good laboratories for addressing this open question of the 
neutrino mass spectrum. 

If non-standard interactions are present then, depending on their coupling strength, the 
picture can change dramatically. In Fig. |2l the shaded region corresponds to the allowed 
values when SUSY FCNC and FDNC interactions are at play. It is obtained by letting the 
A' couplings 8 vary over their entire allowed range - both positive and negative - given in eq. 
([To]) , subject to the further constraints on particular products. 

It is seen that to a significant extent the distinguishability of the two hierarchies is 
obstructed by the interactions unless the number of events is more than about 60. Also, 
the one-to-one correspondence is lost between #13 and the number of events and, at best, a 
lower bound can now be placed on #13 from the observed number. Of course, if the neutrino 
mass hierarchy is known from other experiments, then this lower bound can be strengthened, 
especially for the inverted hierarchy. 

4 The v e — > v T — > r — > [i route is suppressed by phase space for r production and the branching ratio for 
the decay. It contributes at a few per cent level to the signal. 

backgrounds can be eliminated by imposing directionality cuts. The detector is assumed to be of perfect 
efficiency. 

6 This is due to the very tight constraints from \i — > e transition limits in atoms |31| . 
7 Recall we assume that, but for #13 and the mass hierarchy, the other neutrino mass and mixing parameters 
are known. 

8 In fact, we have chosen the subscript m in the A' couplings in eq. I|14|) to be any one of 1,2, or 3. 
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Figure 2: Number of muon events for normal and inverted mass hierarchies as a function of sin 2 2#i 3 
for a five-year ICAL run. The solid lines correspond to the absence of any non-Standard neutrino 
interaction. The shaded area is covered if the A' couplings are varied over their entire allowed range. 



It is also noteworthy that for some values of A'-couplings there may be more events than 
can be expected from the Standard Model interactions, no matter what the value of #13. 
Thus, observation of more than 161 (5) events for the normal (inverted) hierarchy would be 
a clear signal of new physics. 

V.2 Constraining A' 

If #13 is determined from other experiments then it will be easier to look for non-standard 
signals from this f3 beam experiment. However, even if the precise value remains unknown 
at the time, considering the upper bound on #13 one may tighten the constraints on the 
A' couplings. Fig. El reflects the overall sensitivity of the event rate to the $ interactions 
obtained by letting all RPVSM couplings vary over their entire allowed ranges. In this 
subsection, we want to be more specific and ask how the event rate depends on any chosen 
A' coupling. 

At the outset, it may be worth recalling that the BELLE bound on t —* yU7r° [29 severely 
limits the products A 2ml A' 3ml and A 21m A 31m . Thus, i? 23 can be dropped in the effective 
neutrino mass matrix eq. (jUJ). If only A 2ml and/or A' 21m (A' 3ml and/or A' 31m ) is non-zero, then 
R22 (-R33) alone receives an RPVSM contribution. Both i? 22 and -R33 can be simultaneously 
non-zero if A' 21m and A' 3ml (or A 2ml and A' 31m ) are non-zero at the same time. 

In the light of this, we consider the situation where only of the above ^2 coupling is non- 
zero. In such an event, only one of i? 22 , R33 is non-zero. The dependence of the number 
of events on a non-zero A' 331 or A' 2ml , for a chosen sin 2 2#i 3 , can be seen from Fig. El In 
this figure, we use the fact that if only one of these R couplings is non-zero, it appears in 
the results through |A'|. For the normal hierarchy, the curves for A 2ml , for m = 2,3, are 
terminated at the maximum allowed value of 0.18. Fig. |3]can also be used for A' 321 , A 31m and 
A 21m , bearing in mind their different upper bounds. For the inverted hierarchy, the number 
of events is small for A' 2ml and A 21m and insensitive to the magnitude of the coupling. These 
are not shown. It is seen that for the normal hierarchy there is a good chance to determine 
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Figure 3: The number of events as a function of a coupling |A'|, present singly, for the normal (left 
panel) and inverted (right panel) hierarchies. The thick (thin) lines are for |A 331 | (|A' 2ml |, m =2, 3). The 
chosen sin 2 2#i 3 are indicated next to the curves. 



the ^? couplings from the number of events. In fact, if the number of events is less than 
about 50 there is a disallowed region for |A'|, while for larger numbers there is only an upper 
bound. For the inverted hierarchy, more than about five events will set a lower bound on 
the coupling. 



V.3 Effect of A 

The A couplings which can contribute in eq. have strong existing bounds |14j and 

their contribution to R is rather small in comparison to \f2Gpn e - Among them, the bounds 
A121 < 0.05 and A 32 i < 0.07 for m ; ~ = 100 GeV are relatively less stringent P3j. We show their 
very modest impact in Fig. HJ It is clear from this figure that (a) the A-type couplings cannot 
seriously deter the extraction of 9 13 or the determination of the neutrino mass hierarchy, and 
(b) when 9 is is known in future it will still not be possible to constrain these couplings 
through long baseline experiments. 



VI Conclusions 

R-parity violating supersymmetry is among several extensions of the Standard Model crying 
out for experimental verification. The model has flavour diagonal and flavour changing 
neutral currents which can affect neutrino masses and mixing and can leave their imprints 
in long baseline experiments. This is the focus of this work. 

We consider a /3-beam experiment with the source at CERN and the detector at INO. 
We find that the $ interactions may obstruct a clean extraction of the mixing angle 9 is or 
determination of the mass hierarchy unless the bounds on the A' couplings are tightened. 
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Figure 4: Number of muon events for normal and inverted mass hierarchies as a function of sin 2 2^13 
for a five-year ICAL run. The solid lines correspond to the absence of any non-Standard neutrino 
interaction. The shaded area is covered if the A couplings are varied over their entire allowed range. 

On the other hand, one might be able to see a clean signal of new physics. Here, the 
long baseline comes as a boon over experiments like MINOS which cover shorter distances. 
Two experiments of these contrasting types, taken together, can expose the presence of a 
non-standard interaction like RPVSM. 

There are other non-standard models ^H] where four-fermion neutrino couplings with 
greater strength have been invoked. The signals we consider will be much enhanced in such 
cases. 

Our results are presented for the CP conserving case. As #13 is small, the CP violating 
effect is expected to be suppressed. We have checked this for the Standard Model, where 
the 'magic' nature of the baseline [32] also plays a role. 

Finally, in this paper we have restricted ourselves to a /?-beam neutrino source. Much the 
same could be done for antineutrinos as well; then the signs of all terms in R - see eq. (JHJ) - 
will be reversed. It follows from eq. (|21|) that 6"^ can then be maximal only for the inverted 
hierarchy and as such more events are expected here than in the normal hierarchy. Broadly, 
results similar to the ones presented here with neutrinos can be obtained with antineutrinos 
if normal hierarchy is replaced by inverted hierarchy and vice-versa. 
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